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Can Terrestrial Planets Form in Hot-Jupiter Systems? 



Martyn J. Fogg and Richard P. Nelson. 
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Abstract. Models of terrestrial planet formation in the presence of a mi- 
grating giant planet have challenged the notion that hot-Jupiter systems lack 
terrestrial planets. We briefly review this issue and suggest that hot- Jupiter 
systems should be prime targets for future observational missions designed to 
detect Earth-sized and potentially habitable worlds. 



1. Introduction. 



Since the discovery of the first extrasolar planets ( Wolszczan &: Fr"ani 19921 : Mayor k, Queloj 
19951 ). astronomical techniques and observational baselines have advanced to 
the point where o ver 200 extrasolar planetary systems have been identified 



( Butler et al.l 1200^ ). Most detected exoplanets are in the giant planet mass 



range and it is now clear that our solar system is but one variant within a great 
diversity of planetary system architectures. One of the most surprising discov- 
eries has been of a population of giant planets, the so-called hot-Jupiters, found 
orbiting in a region of extreme insolation very close (r < 0.1 AU) to their central 
stars and well within the radius of the ori ginal nebular snowl ine (r ~ 3 — 5 AU) 
where giant planets are thought to form ( Pollack et al.|[l996l ). Hot-Jupiters are 
not uncommon: they amount to about a quarter of exoplanet discoveries, and are 
thought to provide evidence that protoplanets can migrate ov er large radial dis- 
tances via tidal interactions with t he protoplanetary d isk (e.g. lLin Sz Papaloizou 
19861 : iLin et al.l[l996l : IWardlll997l : iNelson et al"]|200 0il. Since the di sk gas is ob- 



serve d to disperse within the first few Myr of the system's existence (jHaisch et al 
200ll ). giant planets must form and migrate through the inner system within 
this period, which is considerably less than the ^ 10-100 Myr thought to be 



required to complet e terrestrial planet fo rmation ( Chamberi 2001 : Kleine et al 
20021 : lHallidavll2004l : lO'Brien et alll2006l 'l. 



Test particle studies have sh own that terrestr ial pla .nets external to a hot - 
Jupiter would have stable orbits ( Jones et al.ll2005l ). and [Raymond et"an ( 20051 ) 
have shown that they should be able to form, in the presence of a giant planet 
already at ~ 0.1 AU, from any available pre-planetary material with a period 
ratio roughly > 3. However, until recently it has been a common assumption 
that terrestrial planets could not have formed in hot-Jupiter systems due to the 
disruptive effect of the giant planet's migration whi ch is deemed to have cleared 
the inner system of planet-forming material (e.g. lArmitagi 20031 ). prompting 
claims that the observed abundance of hot- Jupiters could be used t o constrain 
the general abundance of habitable plari ets ( Ward &: Brownleil2000l ). and even 
their galactic location ( Lineweaver 200ll ). 
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This picture has been chahenged by the work of two groups who have mod- 
eled terrestrial planet forma tion conc urrently with, a nd followin g, an episode 
of giant planet migration (,Fogg fc Nels on 200l, l2006l . l2007a ,b: Raymond et al. 
20061 : iMandell et alJ[2007l ). Their findings suggest that inner system solids disks 
are not destroyed by the intrusion of a migrating giant planet and that terrestrial 
planet formation can resume in the aftermath and run to completion. 

In this paper, we briefly describe our model of terrestrial planet formation 
and show some typical results. 



2. The Model 



We picture our model systems to be composed of a n inner system solids disk, 
at its oligarchic growth stage ( Kokubo &: Idal flQQSl l. extending between 0.4 - 
4.0 AU, with its initial p a,rameters bein g those of a classic minimum mass solar 
nebula (MMSN) model (|Havashil IToSlI ). scaled up in mass by a factor of 3. 
This disk is embedded in the nebular gas which depletes over time by accreting 
onto the central star. We generate a number of different scenarios by allowing 
planetary growth to proceed in the inner disk for between 0.1 - 1.5 Myr before 
placing a 0.5 Mjup giant planet at 5.0 AU and allowing it to migrate inward. 
We end the migration episode when the giant planet reaches 0.1 AU. Since our 
giant is massive enough to open up a gap in the gas disk, we are specifically 
assuming the operation of type II migration, where the planet is locked into the 
visc ous evolution of the gas and migra.tes iiiward in step with its radial mot ion 
(e.g. iLin fc PaDaloizoulll986l : iLin et al.l[l996l : IWardlll997l : iNelson et al"]l2000 'l. 

We use a modified version of the Mercury 6 hybrid-symplectic integra- 
tor to perform our calculations ( Chamberi Il999l ). and run it as an N + N' 
simulation, where we start with N protoplanets embedded in a swarm of N' 
"super-planetesimals" . These latter objects are tracer particles with masses a 
tenth of the initial masses of protoplanets that act as an idealized ensemble of 
a much larger number of real planetesimals and are capable of e xerting realis- 
tic dynamical friction on larger bodies (e.g. iThommes et al.ll2003l ). The central 
star, giant planet, and protoplanets interact gravitationally and can accrete and 
merge inelastically with all other bodies. Super-planetesimals however are non- 
self-interacting but subject to gas drag equivalent to that experienced by a single 
10 km radius planete simal. Details of these aspects of our model are given in 
Fogg Nelson! (|2005l 1. 



We calculate the evolution of the nebular gas using a 1-D viscous disk model, 
with an alpha viscosity of a = 2 x 10~^, that solves numerically a modified 
viscous gas disk diffusion e quation that includes the tidal torques exerte d by 
an embedded giant planet ( Lin &: Papaloizou 1986 : Takeuchi et al.l [19961 ) and 
have described its implementation in iFogg NelsonI (|2007al ). The gas responds 
over time by draining via viscous accretion onto the central star; opening up 
an annular gap centered on the giant planet's orbit; and forming a partial inner 
cavity due to dissipation of propagating spiral waves excited by the giant planet. 
The back reaction of these effects on the giant planet is resolved as torques which 
self-consistently drive type II migration. 
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Figure 1. Details of a migration scenario through a 1 Myr old inner disk. 
The three panels, from the top down, show the system at 20000, 80 000 and 
154 700 years after the start of giant planet migration (ASM). Eccentricity of 
objects is plotted vs. semi-major axis and is read from the left hand axes. 
Black dots are super-planetesimals, grey or white circles are protoplanets, and 
the large highlighted circle is the giant planet. Gas surface density is read 
from the right hand axes, with the upper black curve being the initial profile 
for a SxMMSN and the lower grey curve being the current evolved profile. 

3. An Example Scenario. 

We illustrate the behavior of one of our scenarios in Fig. [H where the system has 
been previously evolved for 1 Myr before the appearance of the giant planet and 
the start of its migration. Three time slices from within the migration episode 
are shown in the three panels: the lowest panel showing the point at which the 
giant planet halts at 0.1 AU, 154 000 years after the start of migration. 

The results show that the passage of the giant does not sweep the inner 
system clear of planet-forming material. Instead, the giant planet shepherds the 
solids disk inward, compacting it and exciting the orbits of objects captured at 
mean-motion resonances. Much of this excited material eventually experiences 
a close encounter with the giant planet and is expelled into an exterior orbit, 
augmenting a new disk of solid material that progressively builds up in orbits 
external to the final position of the hot-Jupiter. In this particular case, 86% of 
the solids disk survives, with 82% of it residing in the external scattered disk. 
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Mass loss is modest and mostly occurs via accretion onto the giant planet. The 
effect of the giant's passage is therefore not the elimination of the inner system 
disk, but instead a modest dilution and strong excitation of solid material and a 
radial mixing that drives volatile-rich material inward. This outcome is robust 
to considerable variation of model parameters, such as the mass of the giant 
and the timing of its formation and migration. The results of further simulation 
of accretion in this scattered disk show that the initially eccentric orbits of 
protoplanets are rapidly damped and circularized via dynamical friction exerted 
by smaller bodies and possibly via tidal drag exerted by the remaining gas 
( Fogg &: NelsonlboOTbl ). Planetary growth resumes and over the following ~ 10 
- 100 Myr gives rise to a set of w ater-rich terrestrial planets in stable orbit s 
external to the hot- Jupiter (see also lRavmond et al.ll2006l : iMandell et al.ll2007f ). 



4. Conclusions. 

Our models predict that terrestrial planets might be routinely expected in hot- 
Jupiter systems, including within their habitable zones, and may be detectable 
by forthcoming missions such as Kepler, Darwin, SIM PlanetQuest and TPF. 
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